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Mean-Square  Error  Due  to  Gradiometer  Field 

Measuring  Devices 

Charles  P.  Hatsell 


Gradiometers  use  spatial  common  mode  magnetic 
field  rejection  to  reduce  interference  from  distant  sources.  They 
also  introduce  distortion  that  can  be  severe,  rendering  experi¬ 
mental  data  difficult  to  interpret.  Attempts  to  recover  the  mea¬ 
sured  magnetic  field  from  the  gradiometer  output  will  be 
plagued  by  the  nonexistence  of  a  spatial  function  for  deconvo¬ 
lution  (except  for  first-order  gradiometers),  and  by  the  high- 
pass  nature  of  the  spatial  transform  that  emphasizes  high 
spatial  frequency  noise.  Goals  of  a  design  for  a  facility  for  mea¬ 
suring  biomagnetic  fields  should  he  an  effective  shielded  room 
and  a  field  detector  employing  a  first-order  gradiometer. 


I.  Introduction 

Measurement  distortions  introduced  by  com¬ 
monly  used  biomagnetic  field  measuring  devices  are 
addressed  in  this  paper.  Comment  is  also  made  on  the 
difficulties  encountered  attempting  to  infer  the  actual 
magnetic  field  from  the  distorted  estimate  provided  by  the 
instrument.  Effects  of  instrument  noise  are  also  discussed. 
Measurement  of  cxtracoiporeal  magnetic  fields  arising 
from  intracorporeal  ionic  fluxes  (viz.,  biomagnetic  fields) 
requires  solving  two  fundamental  problems:  achieving 
high  gain,  low  noise  amplification  and  reducing  the  effect 
of  interfering  environmental  magnetic  fields  (c.g.,  the 
earth’s  magnetic  field,  fields  associated  with  fixed  wir¬ 
ing).  The  former  problem  is  traditionally  solved  by  em¬ 
ploying  as  an  amplifier  a  superconducting  quantum  inter¬ 
ference  device  (SQUID)  [1]. 

Rejection  of  local,  interfering  fields  (which  may  be  6 
to  10  orders  of  magnitude  larger  than  the  biological  fields 
of  interest)  is  in  part  and  almost  universally  approached 
by  using  a  detector  coil  configured  to  reject  low-order  field 
gradients.  The  rationale  for  spatial  common-mode  rejec¬ 
tion  is  that  the  far  field  of  a  source  tends  to  contain  pre¬ 
dominantly  low-order  spatial  components.  The  most  pop¬ 
ular  coil  configuration  is  a  second-order  gradiometer  with 
coaxial,  equal-area  coils  as  shown  in  Fig.  1.  Practically, 
a  second  order  gradiometer  has  beeome  widely  accepted 
as  a  reasonable  tradeoff  between  interference  rejection 
and  instmment  sensitivity.  There  is  an  extensive  literature 
discussing  various  gradiometer  configurations  that  will  not 
be  repeated  here  [2];  especially  relevant  is  a  discrete  spa- 
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Fig.  1 .  Second-order  gradiometer.  Coils  are  of  equal  area  and  coaxial,  with 
the  center  coil  wound  in  opposite  scn.se  to  the  outer  coils.  A.  distance  be¬ 
tween  adjacent  coils  is  the  gradiometer  baseline. 


tial  filtering  approach  to  gradiometer  design  investigated 
by  Bruno  et  al.  [6]. 

A  critical  gradiometer  dimension  is  its  baseline,  A. 
Small  A  enhances  common-mode  rejection  but  decreases 
sensitivity  and  increases  spatial  distortion  in  the  measured 
or  apparent  field  (i.e.,  the  apparent  field  difters  from  that 
field  that  would  have  been  measured  with  a  zero-order 
gradiometer,  also  called  a  magnetometer).  Large  A  re¬ 
duces  distortion  and  increases  sensitivity  -but  at  the  cost 
of  decreased  common-mode  rejection. 

In  any  case,  the  apparent  field  as  measured  by  the  gra¬ 
diometer  is  not  the  actual  field  at  all,  but  some  mildly  to 
severely  distorted  version.  It  is  this  distortion  and  its  mit¬ 
igation  that  are  addressed  here.  The  error  that  finite  coil 
area  introduces  over  a  point  measurement  has  been  ad¬ 
dressed  elsewhere  [3];  so  to  focus  attention  on  the  subject 
at  hand,  it  will  be  ignored  here. 


II.  The  Apparent  Field 

Suppose  our  second-order  gradiometer  is  aligned  along 
the  z-axis  as  in  Fig.  2.  If  the  field  strength  at  the  lower 
gradiometer  coil  at  z  =  0  is  //,(.r,  y,  z),  then  the  apparent 
field  measured  by  the  second-order  gradiometer  would  be 

GlOe,  y,  z)  =  H,{x,  y,  z)  -  2H,(x,  y,  z  +  A) 

+ /f.(.r,  y,  z  +  2A)  (1) 

and  it  follows  directly  that  for  equal  coil  cross-section,  an 
«th-order  gradiometer  would  give 

G'l  =  E  (")(-!)"'//, (.V,  y,  z  +  /«A).  (2) 

»i=0  \tti/ 
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Fig.  2.  Second-order  gradiometer  aligned  with  z-axis. 

While  the  second-order  gradiometer  is  clearly  most  pop¬ 
ular,  the  following  development  will  be  for  nth-order 
equal  coll  area  gradiometers. 

Consider  an  nth-order  gradiometer  oriented  as  shown 
in  Fig.  3.  The  center  of  the  lower  coil  is  located  at  (jc,  y, 
z),  and  the  direction  cosines  of  the  gradiometer  axis  are 
q,  r,  s  giving 

q-  +  +  s'^  =  I  (3) 

and 

u  -  qa^t  +  roy  +  (4) 

where  u  is  the  unit  vector  of  the  gradiometer  axis  and 
Sy,  Oy,  3;  ure  the  coordinate  axis  unit  vectors.  In  a  man¬ 
ner  similar  to  that  leading  to  (2),  the  apparent  field  mea¬ 
sured  by  an  nth-order,  equal  coil  area  gradiometer  is 

G^x,  y,  z)  =  S  +  mqA,  y 

m  =  0  \m/ 

+  mrA,  z  -F  msA).  (5) 


i 


X 


Fig.  3.  nth-order  gradiometer  aligned  with  an  axis  with  direction  cosines 
q,  r,  and  s. 


and  the  fact  that  =  0  (measurement  uncertainties  are 
ignored), 

•  3S.  (11) 

Now  let  n  be  the  unit  normal  to  the  surface  S,  and  define 
a  differential  volume  dV  enclosed  by  5  and  an  external 
surface  S'  everywhere  e  distant  from  S.  Again  from 
Green’s  theorem 

dr  =  '  dS'  -  •  dS. 


(12) 


With  dS  =  n  ds,  and  dS'  =  n'  ds,  an  approximation  for 
(12)  is 


III.  Mean-Square  Measurement  Error 
Throughout,  a  source-free  measurement  space  is  as¬ 
sumed  so  the  magnetic  field  H  is  conservative  and  deter¬ 
mined  by  its  scalar  potential,  <p] 

H  -Vip  (6) 

and 

VV  =  0.  (7) 

A  vector  measurement  error  is  defined  to  be  K,  where 
K  =  H  -G  (8) 

and  the  following  scalar  potentials  are  implicitly  defined; 


G  =  -V<pg 

(9) 

ii 

l 

<1 

(10) 

where  <Pg  is  the  scalar  potential  associated  with  the  appar¬ 
ent  field  G.  If  a  volume  V  of  this  source-free  space  is 
enclosed  by  a  surface  5,  then  from  Green’s  first  identity 


(13) 

which  as  6  0  gives 

V\t  •  Vi/r  *  =  4;  ■  dS.  (14) 

Js  an  Js 

Note  that  the  left-hand  side  of  (14)  is  the  integrated-square 
field  measurement  error  caused  by  measuring  with  a  gra¬ 
diometer  rather  than  a  magnetometer. 

IV.  An  Example 

The  familiar  case  of  a  current  dipole  embedded  a  units 
beneath  a  semi-infinite  volume  conductor  will  be  ad¬ 
dressed.  The  dipole  is  taken  as  y-directed  and  the  surface 
normal  is  n  =  a^.  For  this  conductor  geometry  (14)  be¬ 
comes 

Icortoo 

1  dx  dy  =  \  \  -:rdxdy. 

—  00  J— 00  0Z  V— 00  J— 00  Z 


(15) 
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It  may  be  verified  directly  [4]  that  for  this  geometry  the 
solution  to  Laplace’s  equation  is 


E{a,  /3)  exp  {i2itax  + 


•  exp  [-27r(o!^  +  da  dl3 


(16) 


where  E{a,  |3)  is  the  two-dimensional  Fourier  transform 
of  ^  over  the  ;cy  plane  for  fixed  z  and  is  determined  from 
boundary  conditions;  a  and  |8  are  variables  of  integration. 
Applying  Parseval’s  theorem  to  (15),  the  square  measure¬ 
ment  error  M  becomes 


M  =  J  J  (a^  +  /3^)|£(a,  /3)p  da  d^  (17) 

where 

|£p  =  |Fpll-7p  (18) 

and,  from  Cuffin  and  Cohen  [5], 

F(a,  0)  =  iPya  exp  [-2ira(Q:^  -1-  /3‘)'^V 

[4Tt(a^  -h  f)]  (19) 

on  the  z  =  0  plane;  Py  is  the  current  dipole  moment.  Note 
that  (2)  expresses  the  gradiometer  output  as  a  linear  com¬ 
bination  of  H.  shifted  in  z;  therefore,  in  the  transform  do¬ 
main  there  will  be  some  function  T(a,  0)  that  relates  the 
transform  of  G",  say  g'J  to  the  transform  of  N,,  say  />.  in 
the  following  manner 

g2(«,  0)  =  r(a,  0)  •  /«,(«,  0),  (20) 

Hence  (2)  may  be  expressed  in  the  transform  domain  as 

g'lia,  0)  =  h,ia,  0)  ‘  E 

»i=o  \m/ 

•  exp  [-27rw(a-  +  /3^)'/^A]  (21) 

and  it  follows  from  the  binomial  theorem  that 

T(a,  0)  =  [l  -  exp  (-27r(a-  -1-  0^)'^^^)]".  (22) 

It  is  useful  to  accomplish  (17)  in  polar  coordinates  by  the 
substitutions  a  =  p  cos  6,  0  =  p  sin  6,  and  da  d0  =  p 
dp  dd,  resulting  in 

Joo  p2x 

.  ). 

•  cos^  1  -  Tp  p  dp  de.  (23) 

This  integral  can  be  evaluated  in  closed  form.  The  result 
when  expressed  as  a  fraction  e  of  the  total  magnetic  field 
energy  is 


e  =  I  -  Sia/Af  E  (  ”  )(-l)7(2fl/A  +  mf 
'"=0  xm/ 

2«  /2«\ 

+  4(a/Af  S  (-l)7(2a/A  +  mf 
«'=o  \m  } 

(24) 


n-3 


Fig.  4.  Root-iniegrated-square-measurcment  error  due  lo  a  gradiometer  of 
order  n  and  bareline  A.  Current  source  is  a  dipole  parallel  to  and  a  units 
beneath  the  planar  surface  of  a  conducting  half-space,  e  is  given  by  (24). 


which  is  seen  to  be  dependent  only  on  the  ratio  of  dipole 
depth  a  to  gradiometer  baseline  A  and  gradiometer  order 
n.  Fig.  4  shows  a  plot  of  this  result  for  several  represen¬ 
tative  parameters. 

V.  Magnetic  Field  Recovery 

The  previous  development  suggests  an  inverse  trans¬ 
form  technique  may  be  useful  for  recovery  of  the  mag¬ 
netic  field  spatial  distribution  from  the  gradiometer  output 
spatial  distribution.  This  technique  would  involve  com¬ 
puting  the  spatial  Fourier  transform  of  the  gradiometer 
output  data  and  then  inverting  the  product  of  this  trans¬ 
form  with  \/T(a,  0)  and  an  appropriate  window  func¬ 
tion.  Equivalently,  deconvolution  could  be  performed  in 
the  spatial  domain  by  convolving  the  gradiometer  data 
with  the  inverse  transform  of  l/7(a,  0),  applying  an  ap¬ 
propriate  spatial  window.  Unfortunately,  several  prob¬ 
lems  would  arise  in  such  an  attempt. 

Spatial  deconvolution  without  windowing  is  not  possi¬ 
ble  because  l/Tia,  0)  is  not  Fourier  invertable,  and  even 
with  low-pass  windowing  is  noninvertable  for  gradiome¬ 
ter  order  greater  than  one  due  to  the  singularity  at  the  or¬ 
igin  in  transform  space.  Invertability  would  require  a  win¬ 
dow  that  is  zero  in  the  vicinity  of  the  transform  space 
origin;  this  could  produce  acceptable  results  were  it  not 
for  noise  contributed  by  the  instrument  proximal  to  the 
gradiometer.  This  temporal  noise  process,  which  is  as¬ 
sumed  zero-mean  with  rms  value  «o.  produces  a  spatial 
white  noise  process  of  zero-mean  with  mean-square  value 
%  given  by 

-  («o)^’r  I  \lV(p)/T{p)f  p  dp  (25) 

where  p  =  Va*  +  0^,  and  lF(p)  is  a  windowing  function 
with  properties  as  previously  discussed. 

As  an  example,  the  dipole  model  previously  discussed 
will  be  used.  Let  W(p)  be  a  Hanning  window  modified  as 
follows: 

f  0,  P  <  p„ 

jl,  p,  <  p  <  p,„, 

P  >  P,m 

(26) 


mp)  =  i 


1  +  cos 


vPm 


0, 
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where  is  chosen  to  encompass  a  fraction  3^  of  the  en¬ 
ergy  in  the  current  dipole,  and  =  p„/100.  Using  (19), 
it  is  easy  to  show  that 


baseline=2cm,  dipole  depth=1cm 
peak  field=.5pTjToise=.2pT,  10  averages 


mean  field  true  ■■ —  o'adiometer 

estimate  field  outprtt 


_  _1_ 
"  4ira 


In 


1  +  47ra 
1  -  if 


(27) 


where  the  energy  excluded  for  Pe  >  0  has  been  ignored. 

Since  the  additive  noise  is  zero-mean,  the  expected 
value  of  the  conditional  mean  estimate  is  the  inverse  of 
its  windowed  spatial  transform,  i.e., 

£{^,(x,y)}  =  {{/»:(«. /3)  np) 

•  exp  (i27raA:  -1-  ilv^y)  da  (28) 

where  p  =  Va^  +  |3^,  is  the  two  dimensional  transform 
of  H.,  H.  is  the  estimate  of  H.,  and  E{'}  denotes  expec¬ 
tation.^ 

y)}  can  be  calculated  for  the  dipole  model  by 
using  (19)  and  (26)  in  (28)  and  transforming  to  polar  co¬ 
ordinates  as  was  done  in  (23).  Integrating  first  over  ^  gives 
for  the  z  =  0  plane 


Fig.  5.  True  magnetic  field  ^-component,  its  estimate  using  transform 
techniques,  and  the  output  of  the  gradiomctcr  from  which  the  c.stiniatc  was 
derived.  Measurement  is  along  the  .(-axis  at  :  =  0,  Current  source  is  a 
y-directed  dipole  parallel  to  and  I  cm  beneath  the  planar  (:  =  0)  surface 
of  a  conducting  half-space.  Dipole  moment  is  such  that  the  magnetic  field 
peak  is  0.5  pT, 


irP  V 

£{W,(.r.y)}Uo  =  -rT==^ 

v,v  +  y 

•  p"-^‘''’y,(2irp-'/v^  Hh  y“)p  dp 

(29) 

where  7|  ( • )  is  the  Bessel  function  of  the  first  kind  of  order  p 
one.  This  integral  was  evaluated  numerically  in  the  ex-  ^ 
amples  to  follow;  it  converges  rapidly.  | 

Typical  instruments  used  for  magnetic  field  measure-  | 

ments  have  a  noise  floor  of  about  20  femtotesla  2 

(fT)/'s/Hz,  so  an  instrument  temporal  bandwidth  of  100  * 

Hz  would  produce  an  rms  noise  output  of  about  0.2  pi- 
cotesla  (pT).  Usually,  in  an  evoked  response  paradigm 
some  averaging  is  done  that  for  k  averages  would  reduce 
the  rms  noise  by  1  /■v/^rMaximum  evoked  fields  from  the 
brain  can  be  expected  to' be  about.  0.5  pT. 

Fig.  5  shows  the  z  components  of  the  true  field,  gra- 
diometer  output,  and  expected  field  estimate  as  given  by 
(29),  using  the  modified  Hanning  window  for  a  second- 
order  gradiometer  of  baseline  2'em,  dipole  depth  of  1  cm, 
and  10  averages.  Fig.  6  is  similar  but  assumes  a  dipole 
depth  of  2  cm,. In  each  case  the  temporal  rms  noise  is  0.2 
pT,  (F  =  0.95,  and  the'dipole  moment  has  been  adjusted 
to  produce  a  true  field  peak -of  0:5  pT.  Note  that  even 
though  the  estimate  mean  appears  good,  especially  for  the 
deeper  source,'  ‘the'  rms  'eiror  is  quite  large,  even  for  10 
average^.  Certainly  for  single  event  analysis  of  brain 
evoked  fields,  instrument  noise  would  be  prohibitively 
large,  however,  larger  biomagnetic  fields  such  as  those 
from  the  myocardium  might  be  recovered  effectively  using 

these  techniques.  . * 

As  mentioned  in  the  Introduction,  any  practical  appli 
cation  of  this  technique  would  require  accounting  for  gra- 
diem^ter  coil  area,  because  additional  singularities  are  in 


baseline=2cm,  dipole  clepth=2cm 
peak  field=.5pTjToise=.2pT,  10  averages 


— — •  rriMn  fl»ld  .  tru*  aradionwttr 

Mtimat*  fidld  output 


traduced  into  1  /T{a,  13)  placing  further  constraints  on 
selection  of  p,„.  Finally,  the  reader  is  referred  to  a  recent 
paper  by  Bruno  et  al.  [7]  in  which  an  interesting  but 
somewhat  different  approach  to  spatial  filtering  of  gradi¬ 
ometer  output  is  discussed.  Their  approach  will  also  have 
significant  problems  with  instrument  noise;  however,  they 
deferred  an  analysis  of  the  effect. 

VI.  Discussion 

It  is  clear  that  unless  the  gradiometer  baseline  is  at  least 
several  times  the  source  depth  significant  error  will  result 
from  any  attempt  at  detailed  analysis  of  the  magnetic  fields 
ansing  from  current  sources  of  interest.  Although  many 
investigators  have  without  doubt  been  impressed  with  the 
effectiveness  of  spatial  common  mode  rejection,  the  data 
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collected  have  been  in  many  cases  severely  distorted.  At¬ 
tempts  to  correct  these  errors  through  decorrelation  using 
the  inverse  transform  of  \/T{a.,  /J)  will  not  be  possible 
because  its  inverse  does  not  exist,  and  because  it  is  a  high- 
pass  function  its  use  in  an  inverse  transform  technique 
will  emphasize  high  spatial  frequency  noise  components. 

It  should  be  noted  that  the  explicit  inversion  technique 
explored  in  this  paper  is  model  independent  i.e.,  if  con¬ 
straints  on  the  magnetic  field  structure  are  available,  a 
model  may  be  useful  in  realizing  additional  improvements 
in  signal-to-noise  ratio.  For  example,  some  evoked  re¬ 
sponse  paradigms  used  in  magnetoencephalography  will 
produce  relatively  concentrated  neural  activity  in  the  brain 
[2].  Such  concentrated  activity  might  be  modeled  as  a  di¬ 
pole  source  and  the  data  used  in  a  parameter  estimate  de¬ 
termining  a  best  fit  (e.g.,  least-square)  to  the  dipole 
model.  Such  “additional”  a  priori  information  may  lead 
to  significant  noise  reduction,  and  in  that  sense  an  explicit 
inversion  technique  that  ignores  this  information  incurs  a 
relative  penalty.  Also,  instruments  with  multiple,  spa¬ 
tially  distributed  gradiometers  can  realize  noise  reduction 
through  spatial  processing  techniques.  In  any  case,  the 
inversion  techniques  discussed  here  can  be  applied  at  an 
appropriate  place  in  the  analysis  to  reduce  the  distortion 
introduced  by  the  gradiometer. 

The  analysis  contained  herein  strongly  suggests  that  any 
laboratory  investigating  biomagnetic  fields  consider  a 
room  that  rejects  interfering  fields  and  an  instrument  em¬ 
ploying  a  low  order  gradiometer.  Certainly,  a  first  order 
device  should  be  a  goal. 
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